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Abstract. We consider the first hitting times of the Bessel processes. We give 
explicit expressions for the distribution functions and for the densities by means of 
the zeros of the Bessel functions. The results extend the classical ones and cover all 
the cases. 



1. Introduction. 

In this article we consider the first hitting time of the Bessel process, which itself 
is an interesting object and is one of the important tools to study several problems 
in probability theory. By general theory of one-dimensional diffusion processes, 
the Laplace transform of the distribution satisfies an eigenvalue problem for the 
generator and it is given by a ratio of the modified Bessel functions. 

Except some special cases it is not easy to invert the Laplace transforms. When 
the index v of the Bessel process is a half integer n-|-l/2, n G N, the Macdonald 
function K^, is of a simple form. In this case, it turns that K^^i/K^, is represented 
by the ratio of polynomials. With the help of the partial fraction decomposition, 
Hamana [8, 9] recently has inverted the Laplace transform and applied the results 
to show the explicit form and the asymptotic behavior of the expected volume of 
the Wiener sausage for the odd dimensional Brownian motion. The method used 
in [9] requires some formulae for the zeros of K,^. 

The purpose of this paper is to show explicit formulae for the distribution func- 
tions and the densities by means of the zeros of the Bessel functions. The results 
extend the classical ones, for example, due to Ciesielski and Taylor [4], and cover 
all the cases. 

In order to prove the results we represent the ratio of the modified Bessel func- 
tions by using contour integrals of functions easier to treat, and invert the Laplace 
transform. Recently Byczkowski et al [2, 3] have given similar but different expres- 
sions for the densities of the first hitting times, and applied the results to some 
study on geometric Brownian motions and hyperbolic Brownian motions. We use 
the same curve for the contour integral. But, we show a decomposition of the 
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Bessel function ratio and use it, which makes our expression simpler. Moreover, 
the relation to the classical results in some special cases is clear. 

This article is organized as follows. In the next Section 2 we give the main results. 
After showing two estimates for some ratios of the Bessel functions in Section 3, 
we prove the main results. Theorems 2.1 and 2.2, in Sections 4 and 5, respectively. 
Section 6 is devoted to the asymptotic behavior of the tail probability of the first 
hitting time, which is obtained as an application of the result. In the final Section 7 
we show an addition formula for some ratio of the Bessel function, which is guessed 
from one of our results and is of independent interest. 

2. The first hitting time of the Bessel processes. 
For 1/ e R the one-dimensional diffusion process with infinitesimal generator 



gy") = _ — I- ' _ = : — x 

2 dx"^ 2x dx 2x'^^~^^ dx \ dx 

is called the Bessel process with index v. If 2 z/ + 2 is a positive integer, the Bessel 
process is identical in law with the radial motion of a (2z^-|-2)-dimensional Brownian 
motion. Hence, 2v + 2is called the dimension of the Bessel process. 

The classification of boundary points gives the following information. The end- 
point cxo is a natural boundary for any G M. For ^ 0, is an entrance and 
not exit boundary. For — l<i/<0, Oisa regular boundary, which is instantly 
reflecting. For z/ ^ —1, is an exit but not entrance boundary. For more details, 
see [10, 14] for example. 

For a, 6 ^ we denote by t^I the first hitting time to h of the Bessel process with 
index v starting at a. By general theory of one-dimensional diffusion processes, 
we can evaluate the Laplace transform of the distribution of r'^l by solving an 
eigenvalue problem. In fact, the function 

X I— 7> E[e '^^'f'] 

is increasing (decreasing) on [0,6) (resp. (6, oo)) and satisfies 

0^""^ = Xu, u{b) = 1. 

-A '^"^ 

The following expressions for E[e~-,b] is weU known (cf. [6, 11]): for A > 0, if 
b> and u > -1, 

if < a ^ 6 and v > —1, 

(2.2) ^[e-A^I:.^].«"^^^(«^). 



b-''I^{bV2X) ' 



if < a ^ 6 and v ^ -1, 



(2.3) ^[e-A^i:.^] = «"'^^-^(«^). 



b-''I-^{bV2X) 
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if a > and < 0, 



(2.4) E[e-^"-o] '^"^^ K,{aV2\)- 



if < 6 ^ a and g M, 

(2.5) E[e-<1^, a-'^KAoVTX) 



Here F is the gamma function and Ii, and K^, denote modified Bessel functions of 
the first and the second kinds of order respectively. Both Ijy and K^, are the 
solutions of the modified Bessel difi^erential equation 

(2.6) ,^_+,__(,^ + .> = 0. 

When 6 > and 2i/ + 2 is a positive integer, Ciesielski and Taylor [4] have already 
shown that, for t > 

t 



1 i 

PC^M <t\-i - V '^'^ 



-e 



where is the Bessel function of the first kind of order and {jv,k}'kLi is the 
increasing sequence of positive zeros of Jj^. Hence the density p^^lit) of Tq'^^ is 
given by 

(2-7) Po.M = 2.r(. + i)6^ g ■ 

We will use the notation p^^\ to denote the probability density function of 
While we need to calculate the inverse Laplace transform of the right hand side of 
(2.1) when 2i/ + 2 is not a positive integer, it is possible to obtain it by the same 
methods as those they used to prove Theorem 1 in [4]. Thus we see that (2.7) is 
also valid when 6 > and z/ > — 1. 

In the case a > and < 0, we can easily check 



by (2.4) and the formula 







Formulae (2.2) and (2.5) are found in [1], p. 398, in the case of > 0. Moreover, 
when < a ^ 6, the following formula is provided: 



k 

We shall give a proof under more general situation 
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Theorem 2.1. Let < a < b. Then, for v > -I, 

(2.8) P(.<;.S,).l_2n"£^|^ 



-e 



and, for v ^ —\, 

2u 



-e 
k) 



It should be noted that the result by Ciesielski and Taylor [4] may be obtained 
by letting a — >■ in (2.8) and using the asymptotic behavior of Ju{z) as 2; — >■ 0. 
This theorem immediately shows that, if < a < 6 and v > —1^ 



Similar asymptotic result in the case where a = and 2v + 2is an integer greater 
than 2 (Brownian case) was used in [4] to show the law of iterated logarithm for 
the total time spent by the Bessel process in (0, 6) as 6 J, 0. 

To give our result for the distribution functions of T^l in the case of < 6 < a, 
we need to recall some facts about the zeros of the Bessel function K^. For z/ G M 
we denote by A^(z^) the number of zeros of K^. It is known that N{v) = |z/| — l/2if 
z/ — 1/2 is an integer and that N{u) is the even number closest to \u\ — 1/2 otherwise. 
We remark that N{iy) = if \u\ < 3/2 and N{iy) ^ 1 if |z/| ^ 3/2. Each zero, if 
exists, lies in the half plain G C ; Iie{z) < 0}, denoted by C~. In this case, we 
write Zi,^i, Zy^i^ . . . , ;2i/,Ar(i/) for the zeros. Since is a solution of (2.6), all zeros of 
Ky are of multiplicity one by the uniqueness of the solution of ordinary differential 
equations. This means that all zeros of K^, are distinct. If — 1/2 is not an integer, 
there are no real zeros. For details, see [15], pp. 511-513. 

Theorem 2.2. Let < 6 < a. For ^ and c> 1, we set 

cos{nfi){I^{cx)Kf^{x) - I^{x)Kf^{cx)} 



{K^{x)Y + 7r2{/^(x)}2 + 27r sin(7r/x)i^^(x)/^(x) ' 
(1) If v = ±1/2, 



v2ns^ 



(2) If\u\ < 3/2 and ±1/2, 



b\'^ a — b (g-b)^ 
e 2s 



0/ Jo v27rs^ 



X 



ds. 
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(3) If V — 1/2 is an integer and v ^ ±1/2, 



a 



1- 

Jo V' 



V27rs3 



e 2s ds 



(4) If v — 1/2 is not an integer and Iv] > 3/2, 



-.6 2s cts 



0^/ Jo V27TS^ 

bV^ KJaz^Jb) a-b (g-b)^ I C"-")^ 



2s 6V5 ds 



& \ ^ a — & (g-b)^ 



V2Trs^ 



L\u\.a/b{-l') 

— ' — ^ e ''v^ ax 



X 



ds. 



This theorem gives the asymptotic behavior of P{j^l > t) for < 6 < a, which 
will be discussed in Section 6. It should be mentioned that, when u = 0, it is 
obtained from (2.5). In fact, we have 



,2 10) Jo A K„{bV2\) 

^ 21og(./^ 

Alog(l/A)*-' + ''*-''^ 

as A ^ 0. The Tauberian theorem immediately yields 
(2.11) pri^)>() = «^+„fA 



log t \ log t 

as t ^ cxo (cf. [5], p. 446). In order to derive the second equality of (2.10), we 
have applied the asymptotic behavior of Kq{x) as a; J, (cf. (5.14)). We can 
directly deduce (2.11) from Theorem 2.2 without the Tauberian theorem, however, 
the calculation is left to the reader. 

In case oi v ^ 0, we can not obtain the convenient formula like (2.10) which 
admits us to apply the Tauberian theorem in a straightforward way. 

Remark 2.3. It is well known (see, e.g., [14], p. 450) that the probability laws of 
the Bessel processes with different indices are absolutely continuous. Hence formula 
(2.9) may be deduced from (2.8) and the results in Theorem 2.2 in the case of < 
may be proven from those of u > 0. Although our proofs work in all the cases, this 
remark gives a good check for the results. 
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3. Some estimates for Bessel functions. 

In this section we give two estimates concerning and which we use in the 
proof of Theorem 2.1. Throughout this section, we assume that /x > — 1 and Cj's 
are positive constants independent of the variable. 

We first recall some facts about the Bessel functions. Let D is the set of points 
2; G C \ {0} with [arg^l < tt and {jy^^k^kLi denotes the increasing sequence of 
positive zeros of J^. It is well known that has no other zeros in D (cf. [12], 
p.l27 and [15], pp.478-484). By virtue of the fact that 

(3.1) Jf^ize'^^n = e'^^^J^iz) 

for m e Z and z E D (see [15], p. 75), we see that each —j/^^k is also a zeros of 
for k^l. Moreover, in [15], p. 77, we find 



(3.2) /^(^) 



e-*'^^/2j^(ze^V2) if -TT < arg^ ^ 7r/2, 

gi37rM/2j^(^^g-i377/2) if TT /2 < arg Z ^ TT. 



From this formula we see that the zeros of 1^ with — tt < arg^; ^ tt are jn,k^'^^^^'^ 
for k^l. 

The following estimates for the Bessel functions of the third kind of order 
denoted by H^p and are useful for our purpose. Let 5 > be given. Then it 

is known that, for — tt + 5 ^ arg z ^2t^ — 5, 

(3.3) = ^/^e^^^-'^^/^-'^/^Hl + Ei{z)} 

^ V TTZ 

and that, — 27r + 5 ^ arg z ^tt — 5, 

(3.4) Hi^\z) = ./Ae-»(---M/2-7r/4)|^ ^ 

^ V TT-Z 

where |i?i(z)| ^ C'l/l-^l and |i?2(-2)| = C2/\z\- These estimates can be shown for 
H > —1/2 by the integral representation of the Bessel functions of the third kind. 
Formulae (3.3) and (3.4) are also valid for /x ^ —1/2, since 

H^-li^) = e'^^'-Hj^Hz), H^^liz) = e-'-'^Hj^\z). 

The details are given in [12], p.l21 and [15], pp.197-199. 

The first lemma concerns the exponential growth of the modified Bessel function 

In{z) as \z\ — )> 00. 

Lemma 3.1. Let c and rj he real numbers with < c < 1 and 77 ^ 1. Then there 
exist constants C3, C4 and Ki{c,r]) such that, if R> Ki{c,r]) and Ji^{R) ^ 0, 



< 1 -2(l-c)r7 2 + C3/R 

= 1 - e-2 - C^/R 



holds for every z with \z\ — R and \ axgz\ ^ tt/2 — ATCsm{r] / R) . 



HITTING TIMES OF BESSEL PROCESSES 



7 



Proof. For simplicity we write S for ATCsm{rj / R) . By definition we have 
(3.5) Mz) - l{Hl'\z) + Hl?\z)} 

for 2; e -D (cf. [15], p. 74). We first consider the case where ^ arg2; ^ 7r/2 — 5. It 
follows from (3.1), (3.2) and (3.5) that 



Since -7r/2 ^ arg(^e-*''/2) ^ -5, we obtain by (3.3) and (3.4) that 

1 ^^M/2^(i)^^g-W2) = — ^e^{l + £;i(ze-*'^/2)}, 

2 ^ \/27r2; 

2 ^ V27r2; 

which yields that 

I^,{z) = [1 + £;i(2e-*"/2) + e-^^+^^/^+Vs)-!! + E^^ze-'"" 1'^)} 

Now let z = Re'^ and ^ ^ ^ 7r/2 - (5. Note that 

cos^ ^ cos^^TT ~ '^^ = sin5 = 
Then it is easy to see that 

(3.6) |g-2z±i(M+l/2)7r| ^ g-2flcose < g-2r, < g-2_ 



Hence, if we write 



(3.7) 



I^{z) = |l + e-2-+^(M+i/2)- + ^4(^)1 , 



then there exists a constant C5 such that |-E'4(-2)| ^ C^/lz]. 

Since /^(j^,fce±^-/2) = for A; ^ 1, we have that Iij,{z) ^ under the assumption 
of this lemma. Moreover, by virtue of (3.6), we obtain 



1 + g-2^+i(M+l/2)7r ^ 



> 1 



If > 6*5/(1 — e ^), then the right hand side is positive. Similarly to (3.6), it 
follows that 

I -(l-c)z| < g-(l-c)r? u-2c2;+i(M+l/2)7r| < g-2_ 



8 



Y. HAMANA AND H. MATSUMOTO 



Therefore we deduce from (3.7) that 



(3.8) 



< 



1 



=-(l-c)^ 



|1 + e-2^^- 


i(M+l/2)7r 


+ E^{cz)\ 


|1 + e-2^- 




+ E4{Z)\ 



< 1 ^_(i_e)^ l + e-^ + C5/ci? 



1 - e-2 - C5/R 



for > C5/(l - e-2). 

We next consider the case where —7r/2 + 5 ^ arg^ ^ 0. By virtue of (3.2), (3.3), 
(3.4) and (3.5), we have 

1 



"1 + £;2(^e^'^/^) + e-2^-^('^+i/2)'^{l + Eiize'"'/'')} 



V2 



TTZ 



If we write 



then we have |-E'5(-2:)| ^ C's/l-^l similarly to (3.7). Hence it follows from (3.6) that 
it has the same estimate as (3.8) This completes the proof of this lemma. □ 



Next we show an estimate for the ratio of J^. 



Lemma 3.2. Let c and i] be real numbers with < c < 1 and i] ^ 1. Then there 
exist constants Cq, C7 and K,2{c,r]) such that, if R> K2{c,r]) and Jp(-R) ^ 0, 



< 



2e^"i + 2 + C\JR 
c(2 - Ct/R) 



holds for every z = Re^^ with \6\ ^ Arcsin(?7/i?) and sin(2i?cos0 — nfi) > 0. 

Proof. For simplicity we write 5 for Arcsin(?7/i?). Let z = i2e*^ and |^| ^ S. Since 
< 5 < 7r/6 for i? > 2r], we have that 



(3.9) 



sin^l < sin5 = — . 
' - R 



It follows from |e±»(^-'^M/2-7r/4)| ^ gT-Rsine ^Yia^ 



(3.10) 



-V < |p±»(2-'rM/2-7r/4)| < rj 



for R > 2r]. Therefore we deduce from (3.3), (3.4) and (3.5) that 

(3.11) JJz) = /e*(---M/2-7r/4) ^ g-^(^-^^/2-^/4) ^ ^^^^O ^ 

y/27rz ^ J 
where |i?6(^)| ^ Cs/\z\. Hence we obtain from (3.10) 



1 



(3.12) 



2itsin6> _|_ g-2iisin6l 



+ 2 cos(2if;cos6' - TT/x - 7r/2) + Er{z)}, 
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where |-E'7(2;)| ^ Cg/I^;!- Combining with (3.9), we get 

Finally we show a lower bound of \J^{z)\'^. Note that 

cos(2i?cos^ — TTji — 7r/2) = sin(2i?cos^ — TTji) > 0. 

Using the elementary inequality x + y ^ 2y^xy for x,y ^0, we can obtain 

g2flsine^g-2flsine^2cos(2i?cos^-7r//-7r/2) > 2. 

Hence the right hand side of (3.12) is not larger than 

^-\E7{z)\ > 2 - Cq/R 
2ttR - 2ttR 

If i? > 2/6*8, the right hand side is positive. □ 

4. The first hitting time in the case of < a < 6. 

Prom now on, for a suitable function /, the notation £[/] implies the Laplace 
transform of / and the inverse Laplace transform of / is denoted by All 
formulae concerning Laplace and inverse Laplace transforms can be found in [13]. 
This section is devoted to a proof of Theorem 2.1. For t > and e R let 

A standard formula shows that 



for A > 0. For simplicity we put G^^l{t) = F^^^{2hH). Then we have 



Hence it follows from (2.2) and (2.3) that 

f 1 U{a^^\) 

(4.1) >C[Ga(A) = 



1 I-^{a^) 



iiv> -1, 
if < -1, 



where a = a/b. Note that < a < 1. In order to calculate G^^l{t), we need to 
consider I ^iax) / 1 ^{x) for x > and > —1. 



10 



Y. HAMANA AND H. MATSUMOTO 



Let /I > —1. Recall the formula 



(4.2) I,{z) = Yl 



^n\r{n + ii + l) 



for z & D. For < a ^ 1 and ^ e C we define 

(a^)2^ 



/ \ _ fj, l"^ ) 

<fi^,a[^)-a 2^ 2/^+2nn!r(n + /i + l)' 

n=0 ^ ' 



Since (fiij,,a{ze'^^'^) = 'PiJ.,a{z) for any m e Z, we have that </?^,o is a single- valued 
holomorphic function on C and an extension of z~^I^{az)^ which means that <^p,a(-s) 
coincides with z~^I^{az) for any z E D. Recall that all zeros of are on the 
imaginary line. Hence zeros of coincide with those of including multiplicities. 
This yields that the function ipn^a/ on C is single- valued, meromorphic and an 
extension of I ^{az) / 1 ^{z) . 

For c e (0, 1) and w e D with /^(lu) ^ 0, we set 



and 



for i? > 0, where C(i?) = {z E C ; z = i?e*^, — tt < 9 ^ tt}. The singular points of 
are 0, w and j^^kd^^'^^'^ for /c ^ 1 and they are all poles of order 1. In order 
to apply the residue theorem to S{R), we need to consider the zeros of (fin,! inside 
C{R) and to take R such that a zero of is not on C{R), but it is sufficient to 
put 

Rn = + 

because of the following. It is known (cf. [15], p. 506) that the large zeros of are 
given by the asymptotic expansion 

n+-f^--\ V_i (V _ i)(28//2 -31) 



2"^ 4/ 8(n + /x/2 - l/4)7r 384{(n + ///2 - l/4)7r}3 
In particular, for any e G (0,7r/4) there exists an integer ni ^ 1 such that, for 



(4.3) 

which yields 



(4.4) 



j/i,n -I- £ < -Rn < J/u,n+l — £• 
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Taking the two lemmas in the previous section into account, we take n so large 
that (4.4) and 



Rn > max< \w\,Ki(c, 1), K2(c, 1), — 

hold. Then we deduce from the residue theorem 
E{Rn) = Res(0; /-J + Res{w; /-J 

n n 

+ 5]Res(j^,fee-/2; /-J + Res(j^,fee-^-/2; f-J, 

fc=l k=l 

where Res(f ; /) is the residue of a function / at a pole v. 
It is easy to see that 



^M,i(0) 



Res(0;/;j = = -e 



and that 

Moreover we have 



for each k > 1. Recall 



for z & D. Then, since 



dz 



we obtain 



Resr7 .e^-/2- r ) - - ^/^(cJA^.fcg^"^') 



From (3.1) and (3.2) we deduce 



^M+iOM,fce*''Z2) J,^+i{3^,k) 



and conclude 



(4.6) Ke.O.,e-V»J^--^-|^^ 
Similarly to (4.6), we have that 

(4.7) Res(j^,,e-^'^Z2.j^»j_ ^ Mcj,,k) 
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Therefore, summing up the right hand sides of (4.6) and (4.7), we obtain 



We next prove that E{Rn) tends to as n — )■ oo. We have 
which is equal to the summation of the foUowing three integrals: 



w 



27rX/2 Rnc^'^ - w <p^,i{Rne^^y 



Recall ^^,a{ze^^'^) — V'/i,a(^) for < a ^ 1 and m Then we have 



w f 



'^/^ 1 <^^,e(i?ne^')^^_ 



RnC^^ + W<fii,,i{Rne^^) 

Similarly to Si{Rn), we can show 



= _^r 1 <fi^4Rne^') ^^ 



r/2 

Hence we get 



-/^ 1 UcRne'') 



27r y_,; 



,/2 Rne'^+w I^{Rne'^) 
and 

^^^^ " TT y_,/2 i?^ei2. -^2 /^(i?„e^e) 

Let 

^„ = Arcsi„^. 



de 
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Then < Sn < 7r/2. It follows that S{Rn) is equal to the sum of the following 
three integrals: 



^2{Rn) 



W 



2 r/^-'- 1 I^icR^e 



w 



2 /■7r/2 



1 I^{cRr,e'^) 



w 



2 r-n/2+5r. 1 I^{cR^e'^ 



de, 



7T y_,/2 i?2e^2e_^2 /^(i?„e^e) 



By (4.4), we have J^{Rn) 7^ 0. Applying Lemma 3.1 for 77 = 1 and R = Rm we 
obtain 



\Si{Rr, 



< 



\w\ 



\ 1 ^-2(l-c) 2 + C^/Rr 



TT i?2 — \w 



1 - e-2 - Q/i?, 



(7r-25„). 



This yields that Ei{Rn) tends to as n — >■ 00. 
About S2{Rn)i it holds that 



\^2{Rn)\ 



< 



\W\ 



1 



TT Rl — \w\ 



7I-/2 



/2-5„ 



I^,{cRne 



ie-\ 



Iii{Rne 



de 



\w\ 



7T i?2 — \w\ 







By (3.1) and (3.2), we have 



/^(ci?ne^^e^"/2) 



J^{cR^e 



MRne^') 



de. 



/^(i?ne^^e-/2) 

In order to apply Lemma 3.2 for 77 = 1 and R = Rn, we need to see that 
(4.9) sin(2i?„ cos ^ - tt/x) > 

for 1^1 ^ 5n- It is obvious that 



2Rn COS 6 — TTfl ^ 2Rn — TT/J, = 2n7r + -TT. 

Since 1^1 ^ 5n and sin^n = 1/Rm it is easy to show 



2Rn[ 1-;^) -7TII^2Rr,Jl - ^ -7TllS2RnCOSe-7riI. 



Since i?n > S/tt, we have 



2n7r + ^TT ^ 2i?„ cos 9 — tt/j, 2mr + ^tt 
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and (4.9). Now, applying Lemma 3.2, we obtain 



< 



Rl 



\w\ 



for some constant C\q. This immediately implies that S2{Rn) tends to as n — >■ oo. 
For E'i{Rn)i since 



\Sz{Rn)\ ^ 



\W\ 



1 



TT Rl 



\w\ 



de, 



we can conclude that H'3(i?„) tends to as n ^ oo in a similar way as to E2{Rn)- 
Accordingly S{Rn) tends to as n — >■ oo and the right hand side of (4.8) converges. 
Therefore we have obtained the following. 

Theorem 4.1. //// > — 1 and < c < 1, then 



(4.10) 



li,k) 



for w & D with In{w) ^ 0. 



Recall that ot = a/b G (0,1) and that all zeros of are pure imaginary for 
H > —1. Theorem 4.1 implies that C[G^^l] can be represented by the infinite sum 
of rational functions. Indeed, it follows from (4.1) and (4.10) that, for A > if 
v > -1, 



1 



1 



(4.n) .[«a(A) = ,-^E-^(,,^y,^^^,^^^) 



and that, for A > if ^ —1, 



(4.12) 



a 



2u 



-E- 



In order to prove Theorem 2.1, we need to see that, if ji > —1, there exists a 
positive constant Cn such that 



(4.13) 



Ji^+i{jij,,k) 
for any k ^ 1. By formula (3.11) we get 




2 

TTX 



COs( X — -TT/i — -TT ] + E^^\x] 



for X > 0. Here \E^^\x)\ ^ C^vi jx for some constant ^ which is independent 
of X. Therefore we obtain 



(4.14) 



1 cos(aj^,„ - 7r^/2 - 7r/4) + E^^^^ {aj 



J^+i{j^,n) ^ cos{(j^,, - 7r(// + l)/2 - 7r/4} + 4'^+')(j^,,) ' 



HITTING TIMES OF BESSEL PROCESSES 



15 



Let e e (0,7r/4). It follows from (4.3) that 

(^jn,n - ^TT// - ^TT^ - (n - l)7r 



< e 



for n ^ ni and that 
(4.15) 



cos|j^,n - ^7r{n + l) - ^7r| 



> 



COS£ 



> 



V2 



for n ^ 77-1. We take n so large that 



C 



^ 1, 



< 



2y/2 



and (4.15) hold. Then we obtain from (4.14) 



< 



4\/2 



a 



which implies (4.13). 

We are ready to complete our proof of Theorem 2.1. It follows from (4.3) and 
(4.13) that 



E 



k=l 



jlj.,kJij.+l {jfi,k) 

for each A > 0. This yields 



e-^*e-^M.fe*(it ^ Cii 



'[jf.A^ + jlk) 



< oo 



^ j fj,,kJ^+l{j ^,k) 



j li,kJfi+lij fi,k) Jo 



Therefore we deduce from (4.11) and (4.12) that, if > —1, 



" y.^^Ji',k'Ju+l\Ju,k) 



'''''' e"^-.'=* 



and that, if ^ —1, 



Recalling G'^''i(t) = F^''^{2bH) for t > 0, we complete our proof of Theorem 2.1. 

Since (4.13) assures the termwise differentiation with respect to t in (2.8) and 
(2.9), we obtain the following expression for the densities of r^^J 

Corollary 4.2. Let < a < b. We have that, if v > —1, 



^ Jv+l{ju,k) 



-e 



and that, if v ^ —1, 



- ^ j-v.kJ-u{a3-u,k/ h) 



-e 
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5. The first hitting time in the case of < 6 < a. 

This section devoted to a proof of Theorem 2.2. We again use the same notation 
F^''^ and G^^l{t) as those in the previous section. Set a = a/b > 1. It follows from 
(2.5) that, for A > 



Since Ki, = K^i, for u ^0, it is sufficient to consider the case where u ^0. 

Let 1/ ^ and c > 1. We first assume that z/ — 1/2 is an integer. In this case, 
there is a suitable polynomial i/ji, of order — 1/2 on C such that i/j^iO) ^ and 

(5.1) z'-K^iz) = ^e-^M^) 

(cf. [12, 15]). For example, 

tpi/2{z) = l, 'lps/2{z) = 1 + z, 'ip5/2{z) = 3 + 3z + z'^. 

The function z'^K,^{z) is extended to a entire function and all zeros of i/ji, are the 
same as those of K,^. For 2; e C let 

Then i/j^^c is a single-valued meromorphic function on C and it holds that 

K^{cz) 

= -KM 

for z E D. Therefore, if ;2 G C is not a zero of K^, we have 

K^{cv) 

which implies that K^{cx) / K^{x) can be determined uniquely for a; < if a; is not 
a zero of Ki,. 

Recall our notation z^^i, . . . , -2,y,Ar(;y) for the zeros of K^. Let w is a, point in D 
with Kjjiw) ^ 0. We take R so large that w and all zeros of K^, are inside C{R), a 
circle whose center is the origin and radius R. 

We set 

(5.2) 0{R) = ^ [ 9:;:,Mdz, 

27rz Jc{R) 

where 

9u,c{Z) - 7 X 

z{z — w) 

for z e C The singular points of ^r^^ are 0, w and zeros of K^,, which are all poles 
of order 1. The residue theorem yields that, if N{i') = 0, 

0(i?) = Res(O;^:f,) + Res(w;;<J 
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and that, if 7V(i/) ^ 1, 



N{u) 



G{R) = Res(0; g:^J + Res{w- g^!^,) + Res{z,,r, g^^J. 
By definition of the function ^r^^, we have 



Res(0;^:;;j = -V'.,e(0) = -- 



and 



Res(«;;y-J = e('^-i)-V.,cH = e 



- p(c-l)^ 



. Ki,{cw) 



If N{i') ^ 1, the residue of g'^^^ at Zi,j is equal to 



hm 



for 1 ^ j ^ ^{^)- Since Lemma 3.1 in [9] gives that, if — 0; 

V'i^+i(2;o) = -2;oV'i^(2;o), 

we obtain 

W'll)^{cZ^^j) 



(5.3) 



Res (2;^ J J = - — 



If 2; e D, it follows from (5.1) that 

Ky{cz) ze-^^-^^^i)^{cz) 



(5.4) 



Then Kyicz) j Ki,j^\{z) is extended to a meromorphic function on C. This implies 
that Ky{cx)lKvj^x{x) can be determined uniquely for a; < with Ki,j^\{x) ^ 0. 
Prom (5.3) and (5.4) we deduce 



Res (2^ J J 



and 



e(R) = --+e 



we 



(c-l)z. 



0{R) tends to as — >■ 00 since g^^ci^) = 0{\z\~'^). Hence we obtain 

/ g-(c-l)w 



(5.5) 



Ki,{cw) 



\i V = —, 
2' 



^ ^^j(2i.j - w) K^+i(^^j) ^^^^2 
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in the case where — 1/2 is a non- negative integer. 

We next consider the case where v — 1/2 is not an integer and look for a nice 
expression for Ki,{cvS) j KviyS) Uke (5.5). If v is not an integer, it is weU-known (cf. 
[15], p.80) that 

for z & D and m e Z. When u is an integer, we also have 

K^ize''^'') = e-^'^'^'K^iz) - z7rm(-l)("^-^)^/^(z) 
for 2; e D and m e Z, which is easily seen from 

lim K^{z) = Kn{z) 

for each integer n. Especially, for 2; G -D, we have 

(5.6) K^ize''') = e-'^'^'K^iz) - m/,(z), 

(5.7) K^ize-'"") = e'^'^'K^iz) + ittUz). 

It follows from these identities that the function Kj^[cz) / Ki,[z) can not be extended 
to a meromorphic function on C. For 2 e D let 

^''^^ ^ " ^(^ - «;) K^{z)- 

In order to give a formula for Kj^{cw) / K^{w) like (5.5), we consider the integral of 
^ on a suitable contour. However we can not adopt a circle as the contour like 

(5.2) since k^^^ can not extend to a meromorphic function on C. 
Let e and R be positive numbers with 2e < R. We set 

Or = Arcsin ^ . 
" R 

As a contour, we take the curve 7 defined by 

70 : ^ = Re'^, -Tr + eR^e^TT-eR, 

71 : z = X + ie, —R cos 9r ^ x '^0 

72 : ^ = £e^^ -n/2 ^9^ n/2, 

73 : z = X — ie, —R cos 9r ^ x ^ 
7 = 7o + 7i - 72 - 73- 

We take R so large and s so small that w and all zeros of K,^ are inside 7. Then, 
setting 
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for ^ £ ^ 3, we have 

n{R, £) = ilo + ill - i72 - ils- 
The residue theorem yields that, if N{i') = 0, 

(5.8) n{R, e) = Res{w; /t^^ J 
and that, if 7V(zy) ^ 1, 

Niu) 

(5.9) n{R, e) = Res(w;; /i^^.c) + Yl K,c)- 



It is obvious that 

Res{w;h-J = e ^^^^^ 
When N{i') ^ 1, by using the formula 

(cf. [15], p. 29), we can show 



Res(z,,j;h'^J 



It follows from (5.8) and (5.9) that, if N{v) = 0, 



i7(i?,e) = e(-^)-^^ 



and that, if iV(iy) ^ 1, 



n(R s) = c^c-Dw ^Acw) _ we^^ K,{cz,^j) 

In order to consider asymptotic behavior of n{R,s) as i? — > oo, we recall the 
asymptotic behavior of K^{z) as \z\ — >■ oo from [15], p. 202: if | arg2;| < 37r/2, 

KM = \lf/-'{^ + 0{\z\-')}. 



By combining the formula 



with (3.3) and (3.4), we can make sure that, for | argzj ^ 37r/2 — 5 
(5.10) K,{z) = y^e--{l + E,{ze'^/^)}. 
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Here 5 > is arbitrary given. By virtue of (5.10), we get 

if z = Re'^ e D and 





\w 


1 1 + Ci/c\ 


\z\ 


Vc\z\ 




z — w\ 




z\ 



which tends to as i? — )■ oo. 
For the integral i7i, we have 



ill = — : / 7 r-T- : r — tt—, — — "^C 

2771 J-Rcos 0n + [X + te - W) [x + IE) 

^ .RcosBn ^^c-l){-x+ie) KJc(-X + ie)) , 

ax. 



iJo 



27Ti Jq {x — is){x — ie + w) Ki,{—x + ie) 

Then, using (5.6) and writing the right hand side by 

^,{z) = e-'^''K,{z)-iTTl,{z), 

we get 

w 

27ri Jq (x — is){x — is + w) ^,y{x — is) 

Hence, letting 7^' be the line in D defined by 

7^^ : z = X — is, ^ X ^ RcosOji, 

it holds that 

ill = -— / — dz. 

2711 J^o z{z + W) ^i,[Z) 

Here we define three paths as follows: 

: z = se'^, -7r/2^^^0, 
7^ : z = X, s ^ X ^ R, 
ll : z = i^e*^ -du^e^ 0. 

Since w is inside 7, we have that | Im(t(;)| > £ if Re(t(;) < 0. Recall that there is no 
zero of Kv on the real axis. Then we may apply the Cauchy integral theorem for 
the integral on the contour consisting of 7°, 7^, 7iand 7^ to obtain 



where 



ni = nl + nl-nf, 



nl = — . J d9, 

2m Jg x{x + w) 
1 27rJ_o, Re^^+w URe^') 
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nf tends to as oo. In fact, noting that ^u{xe^^) = K^,{xe^^^+''^) holds for 

a; > 0, we obtain from (5.10) 

^y{Re^) .fc l + £;i(i?e^(^+3-/2)) 

for 1^1 < 7r/6, which yields 



(5.11) 



for large R and a positive constant C13 independent of R and 6*. Since < < 7r/6, 
we see i7f as — )■ 00. 

Furthermore (5.11) shows that the function e~^'^~^^^ ^vic-x) j ^i,{x) is bounded on 
[e, 00) and that Tl\ converges as i? — > 00. Therefore it holds that 

w e-(c-l)ee*'' e (cEe^^\ W g-(c-l)x c ( \ 

R^oo 2-K J_^/2 ee^^ + w Cj.(ee*^) 27^^ a;(a; + ty) ^^.(a;) 

In the same way, we can show that 

fl^oo^ ^ 27ryo ee^^ + iu Ci^(£e*^) 27rz a;(a; + tu) Ci/(a;) 
where 

ioT zeD (cf. (5.7)). Note that 

1 {C,j,{cx) ^iy(cx)1 _ cos{tip){I^,{cx)Kj,{x) - I^,{x)K^,{cx)} 

and recall that the right hand side is Ly^dx). Then we get 

^ .^2 g-(c-l)ee^^ C^(cge^^) ^^ 

_ r°° ^;e-(^-i)^L.,e(^) ^^ 
a;(a; + -u;) 



We will calculate the limit of each term of (5.12) as £ J, 0. 
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Lemma 5.1. Let c > 0, v ^ and \6\ < tt. We have that 

iim —r— — = iim — — — = lim -— — 777— = — 
£4-0 K^{ee''^) eio ^y{ee^'^) sio Ci/(£e ) 

Proof. It is known that 



log0){l + o(l)} 



{1 + 0(1)} ifi.>0 



(5.13) K^iz) 

I 1 

~ 2 

as 1^1 ^> in See [12], p.lll and [15], p.512. Then, it foUows from (5.13) that 

( log(2/rr) - 1 + o(l) .^^^ 
K^(c£e*^) _ I log(2/£) -id 1 + 0(1) 

I 1 + 0(1) 

which converges to 1/c^ as £ 4- 0- Prom (4.2) and (5.13) we deduce that I,y{xe^^) 
converges and K^{xe'^^) tends to infinity as a; 4- 0. This yields that 

Ucee^') K^jcee^') + 0{l) 1 m , nU 

as £ I 0. 

We can show 

^i^ = 1(1 + 0(1)1 

in the same fashion. □ 

The first three terms of the right hand side of (5.12) can be calculated easily. 
Indeed, Lemma 5.1 yields that the first and the second terms converge to 1/40^^ 
and that the third term converges to l/^d" . By (4.2) and (5.13), we can easily see 

^{l + o(l)} if^. = 0. 



r^14^ ^'-.^(x) ^ I (logx, 
■ ^ cos(7rz.) ] c'^(l - c-2-)x2- 

i 2^--W(^ + l) ^^ + ^^^^> 

as X 4 0, which has been noted in [3], p. 29. Hence the last term of the right hand 
side of (5.12) converges as £ 4, 0. Therefore we can conclude 

lim lim n{R,e) = - - / , ^^^^^dx. 



sio R^oo Jq x{x + w) 

Since — for ^ 0, we have that Kyj^i{z) — K\^,\j^i{z) if ;2 is a zero of Ky. 
Moreover, we can regard Z-i,^j as Zu,j for 1 ^ j ^ ^i^)- Therefore we have proven 
the following. 
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Theorem 5.2. Let c > 1, e IR and w is a point in D with K^{w) ^ 0. 

(1) If v = ±1/2, we have 

(2) // \u\ < 3/2 and v ^ ±1/2, we have 



K^(cw) 



-(c—l)'w 



,-(c-l)w 



(3) If v — 1/2 is an integer and v ^ ±1/2, 



L 



x{x + w) 



-dx. 



_ -(c-l)w 



E 



(4) If ly — 1/2 is not an integer and > 3/2, 



-{c—l)w 



y 



we ^)^L|^|^c(a;) 
x{x + w) 



dx. 



We are ready to complete our proof of Theorem 2.2. We have 



1 K^{ay/\) 



a 



a = - > 1. 



We need to invert the Laplace transforms of the following functions: 

Pi(A)=ie-("-i)^, 



The results may be well known (cf. [13]), but we deduce them from the formula 



(5.15) 

At first, put 



2V7rt3 Ja-l 

Then we get by (5.15) 



■ dx = ^Y^ e u,v>0. 



eU-(a-l)}e-^de= ^ 



y/wt Ja-1 



f 

Jo 



e~^Ui{t)dt 



5 3 

e~^^~^t~^dt 



^ / -(«-!)} / e-^-^dx 

Ja-1 Uo 

/•oo 

/ U-(«-l)}e-^^de 

Ja-1 



d^ 



i -(a-l)yA 

A 
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and 

Next we put 



£->i](t) = 9i(t). 



(l2{t;z) = 



Then we obtain from (5.15) 



e-^*q2{t;z)dt 



a-l 
oo 



f 

Jo 

^ / e.U-(a-i)} / 



1 

\/7r 
2 

1 

a/a Ja-l 



4t t 2 



d^ 



,z{e-(a-l)}-^/Ae 



-{a-l)VX 



Vx{Vx - z) 



Hence we get 



>C->2](t) = ?2(t; z). 



Now we have shown, for example, for the fourth case where v — 1/2 is not an 
integer and \v\ > 3/2, 



a 



e 4t 



a-l 



1 1 >^ K^[az^^j) 



e 4t 



a-l 



X 



Finally, a simple change of variables from ^ to s given by ^ — {a — b) ^/2t/s gives 
us the formula in Theorem 2.2 (4). The other cases are simpler. 

6. The tail probability of the first hitting time. 

As an application of Theorem 2.2, we show the asymptotic behavior of P{t^I > t) 
as t — )■ oo when < 6 < a. In Section 2 we showed it when = by the Tauberian 
theorem. In [3], it is shown that, if < 0, 

P{ri:i>t) = cj'^ {1 + 0(1)} 

holds for some constant c,y. It should also be noted that, in [16], Yamazato has 
discussed on the tail probability in a general framework, and some Bessel processes 
may be treated. We give an explicit expression for the constant c^. 

To make the statement clear, we define two constants when — 1/2 is an integer. 
Put 



(^) _ (a-b) 



2M 



0■^ = 



2\v\ 
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Moreover we set CTg'^^ = if ji^l = ±1/2 and 

if otherwise. 

Theorem 6.1. Let < 6 < a. 

(1) 7/^^ = 0, 



log t \ log t 

(2) If v > and v — 1/2 is an integer 



» . ,^ _ , r ''^ ^ f>>\ /2(-l/2r-'/^f J.) ^ /aV Ml 1 



ay ' \aj V TT (zy-1/2)! 1^"^ ' V^>y ^ j 



(3) If v < and v — 1/2 is an integer 



(4) If v > and v — 1/2 is not an integer, 

(5) If ly < and v — 1/2 is not an integer, 



Remark 6.2. It seems that (4) and (5) also hold when v — 1/2. But we do not 
pursue the identities. 

Before proving this theorem, we give two lemmas. The following one is the 
immediate consequence of Lemma 4.3 given in [3]. We let m{v) be the greatest 
integer which is not larger than \v\ — 1/2. 

Lemma 6.3. We assume that \v\ ^ 1/2. If v — 1/2 is an integer, we have that, 
for any ^ m ^ m{v) — 1 

(6.1) lim t"^+^/2p(Ti"2 > = 0. 

If V — 1/2 is not an integer, we have (6.1) for any ^ m ^ m{iy). 
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For t > 0, IV 7^ and z E C~ set 



vStts^ 

* a — b (a-b)^ 
e 2s 



V^TTS^ 



X 



ds. 



Theorem 2.2 implies that P{t^I > t) is represented by a hnear combination of ^i^s. 
Changing variables from s to tt given by {a — b)/ y/s = u, we have 

/2 r{a-b)/Vt 2 



flit) = 1 



e '2 du 



e 2 du. 



71" J{a-b)/Vi 



^ J{a-b)/Vi 



e ^ ^ b du, 



where 



J(a-b)/Vt UO 



' ' ' e b dx 



X 



du 



and L^\iy\ = L|,y|^a/;,(a:)/ cos(7rz/). It is obvious that Lj^^| is positive on 

(0, oo) since I^, and K^, is increasing and decreasing on (0, oo), respectively. For an 
integer m with ^ rn ^ m(zv), we set 

(a - 6)2^+1 



/Si(m) 



2m + 1 ' 



1 f a — b 



P^^\m)^{2m)\b^-^^'Y.k\\ b 

fc=0 



e b dx. 



^ 



2m-fc+2 



2m+l 



^(g-b) 1 J z{a — b) 



e b 



E 



fc=0 



ix) = n\i -e 
a 



^(a/b-3)x 



{1 + 0(1)} 



where 

(3{m; z) = -(2m)! I 
It follows from (3.7) and (5.10) that 

as a; — >■ OO. Moreover, (5.14) yields that L^^^ a/b(^)/^^"^~'^^^ asymptotically equal 
to a constant multiple of for small x. Since 

2m - /c + 2 - 2|iy| ^ 2m{u) + 2 - 2|zy| < 1, 
we have that the following improper integral converges: 

—hr — , , „ e dx = cosy'Kv) 

r.2m— fe+2 ^ ' 



f 

Jo 



X^ 



e b dx. 



X 



2m-fe+2 



HITTING TIMES OF BESSEL PROCESSES 



27 



Lemma 6.4. If \iy\ ^ 1/2, 

-m{i/) 



^(f^ (-l/2)-^i(m) 1 , ^/ 1 \ 

m=0 ^ ^ 

(6.3) tf^2(t, ^) - V + ) 

m=0 ^ ^ 

If 0<\iy\< 1/2, 



/f > 1/2 and — 1/2 is not an integer, 

[2-^ {-l/2rpt\m) 1 



(6.5) (h- \ ' \ I a 



m=0 

2 \ I'^i r \ 1^ 

+ I ^ 



+ 



2 y \V^/ va/ J r(i + |iy|)ti^i 
1 



Proof. For a; ^ let 

m=0 ^ ^ 

Note that 
Hence we have 

AT '5^) 1 / f{a-b)/Vi fj i'{a-b)/Vt 



which imphes (6.2). Similarly, by the formula 



poo ^ 

(6.7) / a^^e-'^^dx = e''^'' V 



n 



A;! n^-k+^ 



for /3 > and Re(/x) > (cf. [7], p. 340), we immediately get (6.3). 
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For 1 1^1 > 1/2 and v — 1/2 is not an integer we have 



m=0 



m\\ 2 



u 



2m 



{a-b)/Vt 



X 



'{a-b)/Vt 



f 

Jo 



i/b 
X 



e b dx 



e dx 



du. 



du 



Then the first term of the right hand side is equal to 



m=0 



cos(7rj^)m! pn+i/2 



since 



J{a-b)/Vt UO 



e dx 



X 



du = 



which is obtained by the Fubini theorem and (6.7). We set 



J(a-b)/Vi Uo 



' ' ' e b dx 



X 



l{a-b)/Vt 

Changing variables from x to y given by xy/iu/b = y, we have 



du. 



t-oo r /"OO 



f_by_ 
'^'Wtuj y 



dy 



du, 



where 1a is the indicator function of A. To see the convergence as 
t — > OO, we need to dominate 



(6.8) 



u/ y 



by an integrable function which is independent of t. We have that (6.8) is equal to 



(6.9) 
Since 



{u)\P^''\u) 



Llua/biby/Vt 



u 



{by/Viu)^\'^\ 



_y2H-l^-2\.\^-y_ 



-^H,a/b(^) ^_(a/b-3)x 



x 



2W\ 



is bounded on (0, oo), we have that (6.9) is dominated by a constant multiple of 



(6.10) 

We have that, if a < 36, 



(«)|p('^)(«)|e'^-Vl^|-'«~'l"l 



y (a — 3b) 

e Vtu ^ < 
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and that, if a > 36, 

g Vtu \ g a — o ^ g a — o^ 

for {a — b)/ Let 

K = mm< 1, 



a — b 



and hence (6.10) is bounded by 

(6.11) |p(^)(«)|«-2|'^ly2|^l-ie-''f. 

To see that \P^'^\u)\u~'^\'^\ is integrable on (0, oo), we note that 

\P^''\x)\ ^ Ci4min{l,a;V2}a;2"»(^) 
for some constant C14. Then we get 

|p(^)(a;)|«-2|'^l ^ 

Since 



Ci4w2™M-2|^H-2 if < W ^ 72, 



2m{v) - 2\v\ + 2 > 2( - ^ ) - 2\v\ + 2 > -1, 









-2|zy| 







we see that the function given by (6.11) is integrable on (0, 00) x (0, 00). Applying 
the dominated convergence theorem, the Fubini theorem and (5.14), we have that 
t\''\^^{t;v) tends to 

i"-!^) 22iH-ir(|H)r(|H + i) Jo ' ' io ^ * 

as t ^ 00. Since 

(cf. [12], p.3), (6.12) coincides with 



Changing variables from tt to f given by f = tt /2, we have 

[ P<')(„)„-Mrf„ = ^ j- _^|,-« _ 'g' tip.™}*. 

^ ^ ^ m=0 
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which is equal to 

1 :r(l-M 



21^1 + 1/2 V2 

(cf. [7], p.361). The formula 



2 / \2 / cos(7rz^) 
(cf. [12], p. 3) immediately yields 



2 cos(7rzy) V 2 / {\b J \a J jr{l + \iy\)&\ \& 

and hence, we have (6.5). 

When < Izvj < 1/2, it is enough to consider W^{t;i') directly. We can easily 
deduce (6.4) in the same way as W^{t;i') for lu] > 1/2. The calculation is left to 
the reader. □ 

We are now ready to prove Theorem 6.1. We need only to show Theorem 6.1 in 
the case of 7^ 0. 

For simplicity we set c = 6/ a and 

Wj = ^ — — — 

for 1 ^ J ^ ^{^)- We first consider the case when u — 1/2 is not an integer and 
u <0. In this case, if < —3/2, Theorem 2.2 gives 

Piri:i >t) = Mt) + E ^i'^'^2(t; z.j) + c^Ut; ^). 

Note that m(z/) + 1/2 < < m(z/) + 3/2. It follows from Lemma 6.4 that 

m=0 

By virtue of Lemma 6.3, we have that, for any ^ m ^ m(i^) 

(6.13) A (to) + c"/?^"^ (m) + c"4"^ (m) = 0. 
This immediately yields 

(6.14) P(rS > t) = (A)^c- - c-)^ + 0(e). 
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If -3/2 < < 0, Theorem 2.2 gives 

In the case of —3/2 < < — 1/2, we have by Lemma 6.4 that 

p(4:^ > «) = yflvm + <rpi"\m + (^jV" - --)Ti^) + 

Lemma 6.3 yields 

A(0)+c'^4'^^(0) = 0, 

and hence we have (6.14). In the case of —1/2 < z/ < 0, Lemma 6.4 gives (6.14) 
since ^i{t) is of order l/\/t. We therefore obtain Theorem 6.1 (5). 

We next consider the case when — 1/2 is not an integer and v > Q. Ifz^>3/2, 
it follows from Theorem 2.2 that 

i=i 

which is equal to 

m=0 

V 2 y ^ ^r(i + z/)t^ V^'' / 

We have that (6.13) is equivalent to 

for ^ m ^ m{y). Since ^-jyj is regarded as z^j,j and K-jj^\{z-jj^j) — K^jf-i{ziyj) 
for 1 ^ J ^ ^iy)^ we obtain that I3^\m) = ^\m). Moreover it is trivial that 
P^\'m) = ^2 '^\'m), and hence we have that, for ^ m ^ m{iy) 

Pi{m) + c-''f3i''\m) + c-''f3i''\m) = 0. 
This immediately yields 

(6.15) P(Ti") >t) = l-c^''+(^] \c-^ - c'')-—J—— + o(-). 

In the other cases, (6.15) can be derived in a similar way. Then we finish to prove 
Theorem 6.1 (4). 

We lastly consider the case when u — 1/2 is an integer. Similarly to the case 
when u — 1/2 is not an integer, we can deduce the asymptotic behavior for z^ > 
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from that for v < 0. Hence we shall treat only the case of < 0. If v —1/2, 
Theorem 2.2 and Lemma 6.4 give 

m=0 ^ 

Note that m(^') = —v — 1/2^1. It follows from Lemma 6.3 that 
Since 

if = —1/2, we obtain Theorem 6.1 (3). We finish to prove Theorem 6.1. 

7. An addition formula for the sum of Bessel functions ratio 
Prom the results in Section 4, we have that, if > — 1, 

E\e-^<h = 1 - 46^ f - V V - Maj.,k/b) 



bY^f 1 jjy,k \Jv{a3u,k/b) 



Hence, heuristically, we obtain the following addition formula, which may be of 
independent interest. In fact the idea of proof is the same as that in Section 4. 

Theorem 7.1. Let > — 1 and < c < 1. Then we have 



^-^ j^^,nJ^J,+ l{j^J.,n) 2 



Proof. Letting (p^^c be the holomorphic function considered in Section 4, we put 
for 2; e C with (fin^i{z) 7^ and 



^TTi Jc{R) 
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for R > 0. We take n so large that (4.4) holds. Hence the residue theorem yields 
that 

n n 

S{Rn) = Res(0; l^,c) + ^ Res(j^,fce^'^/2; l^,c) + Res(j^,fce-*"/2; l^,c). 

fc=l k=l 

It is easy to see 

Res(0;/^,c) =c^. 
We can show, in a similar way to (4.5), 

Eesd ^e^^/2•/ ^- -^t^^^^ti^) 

and, in a similar way to (4.6), 



JlJ,,k'Jfi+l \JiJ.,k) 



Therefore it follows that 



k=l ^l^'kJix+l{jii,k) 

We now prove that S{Rn) tends to as n — > oo. We have 



f^^ARne'^) 

which is equal to the sum of the following three integrals: 

^l(-Rn) = 7^ / '^^'' /p" jfl! ^^' 



r/2 ^ix,l{Rne^^)' 

Since (/?^^c(^e^*'^) = (/?^^c(^): the change variable formula yields 

2^i[Rn) - 2^ <^^,i(i?„e^(^-))^^ - 27r io /^(^ne^^) 

and 

2^s[Rn) - 2^ ^^,i(i?„e^(^+-))''' - 27r I^/^ /^(ii^e^^) 

Hence S2{Rn) = S^{Rn) + ^siRn) and 



(7.1) \mn)\ ^ I 





/^(ci^ne^') 


/-7r/2 


/^(^ne^^) 



d9. 
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Let rj ^1 he given and set 



Sn = Arcsin 



Rn 



Note that 5„ tends to as n — >■ oo and that the right hand side of (7.1) is equal to 



(7.2) 



l>n/2 — S„ 




'-n/2+Sr, 


URne^') 



IT 



tt/2 



n/2-6r, 



I^^iRne' 



de 



+ 



/.-7r/2+5„ 



TT 



7r/2 



I^{Rne'') 



de. 



For the first term, by using Lemma 3.1, we get 



/7r/2— 5„ 
-7r/2+5^ 



i0-\ 



URne^') 



de < Ci^e-^^^-"^" 



for some constant C15. From (3.1) and (3.2) it is easy to see for the second term of 
(7.2) 



.7r/2 



URne^') 



de 



I,{cR^e^('+-/^y 



lJRne^iS+-/^)) 



de 



MRne^') 



de 



and from Lemma 3.2 





J^{cR^e'') 




MRne^') 



de S Ci6e^^5^ 



The third term of (7.2) is estimated from above in the same way. 
Hence we get 

hmsup|r(i2„)| ^ Ci5e-2(i-'^)'? 



for an arbitrary ry ^ 1 and, therefore, U{Rn) tends to as n — )■ 00. □ 
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